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Abstract. Broadband femtosecond-laser frequency combs are filtered to spectrographically resolvable 
frequency-mode spacing, and the limitations of using cavities for spectral filtering are considered. Data 
and theory are used to show implications to spectrographic calibration of high-resolution, astronomical 
spectrometers. 
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1 Introduction 

High-resolution astronomical spectroscopy has seen rapid 
advancement in the last decade. Increased light gather- 
ing and instrument stability suggest the possibility of ra- 
dial velocity measurements of astronomical sources with 
lOcm/s precision. This precision would allow observation 
of terrestrial- mass planets, measurement of temporal vari- 
ation of fundamental constants, and direct observation of 
the acceleration of the expansion of the universe. The 
inherent limitations of existing astronomical wavelength 
standards (iodine cells and Th-Ar lamps) inhibit attaining 
the level of precision necessary for these measurements. An 
ideal wavelength standard would provide a high-density 
array of uniformly spaced, constant spectral brightness 
emission lines, the frequency of each tied to fundamental 
constants or to the standard SI second. In addition, the 
precision and long-term stability of this source should ex- 
ceed the baseline of a high-resolution spectrograph. Optical 
frequency combs based on femtosecond lasers address all 
of these requirements |l|2j . 

The current generation of frequency combs produces 
an array of evenly spaced frequency modes spanning hun- 
dreds of nanometers [3,4^ . While numerous types of fem- 
tosecond laser combs exist, the spectral range of astro- 
nomical interest from ~400-2000nm is well covered by two 
types of passively mode-locked lasers that employ mature 
technology [5 and that have been frequency-stabilized. 
The first type is based on Ti-doped sapphire, and the 
second class centers around Er-doped optical fiber tech- 
nology. Tiisapphire femtosecond lasers with intracavity 
spectral broadening centered at 800 nm have a typical 
fractional bandwidth AX/X ^ 0.75 (coverage from ~600- 
1200nm) ^6f7|l8j . Still greater spectral coverage, 400-1500nm, 



can be achieved by spectrally broadening the output of 
the Tiisapphire laser in nonlinear microstructured optical 
fiber [9|10|llj . In the near-infrared, femtosecond Erifiber 
lasers are readily broadened in highly nonlinear fiber to 
provide coverage from 1000-2200 nm [12 13 14 . Despite 
sufficient spectral coverage by these lasers, the small fre- 
quency mode spacing of ^IGHz for Tiisapphire and <250MHz 
for Erifiber prohibits resolution of individual frequency 
modes in typical astronomical spectrographs. 

This paper discusses Fabry-Perot cavity filtering of 
femtosecond-laser frequency combs, ultimately resulting 
in a wide spectrum of resolvable frequency modes with 
line spacing in the range of 10 to 30 GHz. Cavity filtering 
is demonstrated by use of a ~ 1 GHz Tiisapphire laser, 
and the effect of mirror parameters, refractive index of 
the cavity, offset frequencies, mode suppression, and comb 
linewidth are explored. The experimental results are com- 
pared with a theoretical model, and the implications of 
using a comb as a spectrographic calibration standard are 
considered. 



2 Cavity Filtering Theory 

Using a /rep^^l GHz, Tiisapphire comb as a prototypical 
laser, the basic principle of cavity filtering is illustrated 
in Fig.[T] A mode-locked laser produces a comb of lines 
equally spaced in frequency [15], which in current, state- 
of-the-art systems are spaced too tightly in frequency to 
be resolved by standard astronomical spectrographs. As 
suggested by Sizer |16] , a Fabry-Perot cavity may be used 
to spectrally select a set of modes from the comb that 
have a frequency spacing larger than /rep. The length of 
the cavity is adjusted such that each cavity peak passes 
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Fig. 1. A /rep= 1 GHz mode-locked laser is filtered to a 
frequency-mode spacing of mf rep by use of a Fabry-Perot cav- 
ity. The light reflected from the cavity is picked off by a polar- 
izing beam splitter (PBS), passes through an optical, band- 
pass filter (fl), and is the input to a PZT-based servo system 
controlling the length of the cavity. The lower section shows 
cavity amplitude (blue dashed) and phase (green dotted) that 
act on the electric field of the comb. An offset between the 
cavity and the comb at the middle comb mode is defined as 
Aucc = z^fp_mid — z^c_mid- The modes to the right and left of 
the m*^ modes are defined as the right (nur) and left (nni) 
nearest-neighbor modes. 



a comb mode that is separated from the previous cavity 
mode by m comb modes, thereby creating a comb with 
frequency spacing m times that of the input. The greater 
the mirror reflectivity, the higher the suppression of the 
off-resonant modes. 

This straightforward concept of using a cavity to sup- 
press intermediate frequency modes is complicated by hun- 
dreds of nanometers of spectrum necessary for spectral 
calibration |17|18j . While the frequency modes of the laser 
are equally spaced across the entire spectrum, the cavity- 
mode spacing is frequency dependent. Dispersive media 
inside the cavity as well as additional phase contributions 
of mirror coatings cause the frequency mode spacing or 
free spectral range (FSR) of the cavity to vary across the 
spectrum. If the comb and cavity are well-overlapped at 
one spectral region, they walk off one another in other re- 
gions. In order to maximize the bandwidth of light trans- 
mitted through a cavity while maintaining suppression 
of neighboring frequency modes, the mirror reflectivity 
R{uj), mirror phase (/>r(^), refractive index inside the cav- 
ity n(a;), and shift frequency Aucc{^) between the cavity 
and the comb must be considered. The offset frequency 
Aucci^) is defined as the frequency shift between the peak 
of cavity transmission and the nearest comb mode at a 
given frequency. 

In designing a large-bandwidth filter cavity that ex- 
hibits high suppression of intermediate frequencies, the 
most significant terms develop from the mirror coatings 
themselves. For low- group- velocity-dispersion mirrors such 
as those in Fig.[2j there is a trade-off between higher reflec- 
tivity (which increases suppression of neighboring modes) 



and the bandwidth over which the cavity has a uniform 
FSR (where the cavity can be well matched to the comb). 
Even small changes in the FSR are detrimental to the us- 
able spectral bandwidth of the cavity. Like a Vernier scale 
with one metric having inconsistent spacing, the walk-off 
between the comb and the cavity is cumulative, with each 
successive m^^ comb mode slipping further from the cavity 
mode. 

At first inspection of the mirror reflectivity in Fig.[2j 
it appears that the mirrors' coatings have bandwidths of 
a few hundred nanometers. When accounting for mirror 
phase, the bandwidth over which the cavity can filter the 
comb is much less than the region with acceptable re- 
flectivity. Assuming the mirror phase is linear in wave- 
length, the index inside the cavity is unity, and R ^ 1, 
the maximum mirror phase that can be tolerated while 
maintaining lout/hn ^ 0.5 is |^| > ^/{R^^l^^JR. Note 
that ^ refers to the difference in mirror phase across the 
total usable bandwidth beyond the linear phase acquired 
due to normal propagation. In actuality, the mirror phase 
may oscillate between positive and negative values until 
the cumulative phase exceeds that of the above expres- 
sion. Here we can neglect the Gouy phase shift, which has 
negligible frequency dependence. Using the above expres- 
sion and the mirror phase in Fig.|2j a cavity bandwidth of 
~ 100 nm centered at 800 nm should be achievable. 

When the cumulative mirror phase is small, such as 
for specially designed coatings or for metal mirrors, the 
index of refraction of air inside the cavity can become the 
dominant term affecting usable bandwidth. Assuming a 
linear variation of the index as a function of wavelength 
within the cavity (a reasonable approximation for air in 
the visible spectrum) the bandwidth A\ over which comb 
transmission ^ 50 % is limited by the change in refrac- 
tive index {An) as A\ < frei>\/{R - lY / R/ (ttcAti), 
where A is the center wavelength of the cavity. The filter 
bandwidth varies inversely with An] whereas, the depen- 
dence on /rep and m is linear. This somewhat counterintu- 
itive linear relation stems from /^ep and m being inversely 
proportional to cavity length. A refractive index change 
of up to 1 part in 10^ allows an optical bandwidth of the 
cavity of 100 nm. 

The above equations give a terse understanding of the 
cavity limitations. To take properly into account the cavity- 
comb interaction, the system is modeled by considering 
the electric field of all delta-function comb lines, each 
spaced by a frequency /rep, interacting with an air-spaced 
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Fig. 2. Mirror coating data provided by Ramin Lalezari of 
Advanced Thin Films. The phase data include the linear prop- 
agation component that is subtracted for use in calculations. 
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Fabry-Perot cavity of nonabsorbing mirrors. The steady- 
state, planar or tilted- wave cavity equation is used to 
model the effect of the cavity on the amplitude and phase 
of each comb mode [19 : 

1 _ 7^(^)ei[2nMa;L/c+20,(a;)+0Z5] ' V J 

R{(jo) and ^r(^) account for the coefficient of reflection and 
phase upon reflection of the cavity mirrors. Because this 
is a planar derivation, is added empirically and con- 
siders all geometrical phase terms such as the Gouy phase 
(~ 2cos~^[l — L/r] where r is the radius of curvature of 
the cavity mirrors). The index of refraction, n(a;), is taken 
for that of air using the Ciddor equation [20 for typical 
laboratory conditions of 24 °C, 630 Torr, a fractional hu- 
midity of 30%, and 400 ppm of CO2. The cavity length 
L is adjusted such that the filter cavity passes an inte- 
gral number m of comb modes through each cavity FSR. 
Unless otherwise noted, all theoretical curves assume that 
the comb modes overlap the cavity modes; Av^^ = at the 
center of the mirror coatings uOmid- 

Modeling the cavity parameters by use Eq. [l] repro- 
duces the observed spectrum in both the time and fre- 
quency domains. Although cavity filtering is practical for 
increasing the temporal pulse rate |16|21|22] as well as for 
improving performance of microwave generation by opti- 
cal frequency combs [23"'24[, here we focus on the cavity 
effect of thinning comb frequency modes. By taking into 
account the amplitude and phase of the cavity from Eq. [l] 
on each of the incident comb frequency modes, a theoret- 
ical model of the output spectrum is possible. 



3 Experimental setup 

Figure [1] depicts the apparatus used in the experiments. A 
broadband, Ti: Sapphire laser [25] generates a ~lGHz rep- 
etition rate, octave-spanning (550nm-1100nm) comb. The 
frequency of the N^^ cavity mode may be written as /at = 
N frep + fo' The comb has two degrees of freedom: the 
carrier-envelope offset frequency fo and the mode spac- 
ing /rep. The repetition frequency is varied with cavity 
length and is easily accessible and controllable. The offset 
frequency is less readily determined. Taking advantage of 
the full comb spectrum which exceeds an octave of band- 
width, a group of modes at one end of the spectrum may 
be frequency doubled and then compared to comb modes 
on the other side of the spectrum. By adjusting the pump 
laser power with an acousto-optic modulator, the observed 
beat frequency which corresponds to fo may be referenced 
to a known frequency. Both fo and /rep are ultimately ref- 
erenced to a hydrogen maser; the resulting spectrum con- 
sists of ~ 10^ comb frequencies, equally spaced by /rep. 
To more closely match the spectrum of the lasers to the 
reflective center of the cavity, the laser is first filtered by 
multiple reflections from plane mirrors of the same coating 
as the cavity. 

The remaining spectrum is mode-matched into a Fabry- 
Perot filter cavity of variable length. The cavity consists 




750 800 850 900 950 1000 1050 



Wavelength [nm] 

Fig. 3. (a) Experimental (green solid) and simulated (grey 
dashed) transfer function (/out/^in) for the m = 20 comb 
modes. Cavity length L ^ 0.75 cm is locked with filter fl 
centered at 980 nm. (b) Effect of cavity lock wavelength (as 
determined by fl) on transmission through the cavity. Cavity- 
filtered comb spectra in power per mode as a function of wave- 
length are shown for the lock wavelength centered at 820 nm, 
840 nm, 860 nm, 880 nm, 900 nm, and 960 nm. The input spec- 
trum has been vertically scaled by m = 20 to account for the 
suppression of intermediate modes. 

of two spherical mirrors with radii of curvature 50 cm and 
dielectric coatings of peak reflectivity 99.2 % centered at 
910 nm. In order to stabilize the cavity length, the light 
reflected from the cavity is separated by a polarizing beam 
splitter and is detected on a photodiode. A tunable, opti- 
cal bandpass filter (fl in Fig.jl]) selects the portion of the 
optical spectrum used to lock the cavity. The detected 
light creates an error signal, which maximizes the inten- 
sity of the cavity-transmitted modes in the spectral band- 
width of fl. This error signal drives a PZT servo system, 
dithered at 74 kHz, that controls the length of the cav- 
ity. The cavity length may be scanned and then locked 
such that it passes m comb modes per FSR. Because the 
Tiisapphire laser is not a purely TEMqo mode, excitation 
of only the lowest-order cavity mode requires additional 
spatial filtering. The comb teeth that traverse the cavity 
in the lowest-order cavity mode are detected either with 
a fiber-coupled fast photodetector to obtain a RF spec- 
trum or with an optical spectrum analyzer to measure the 
optical spectrum. 



4 Experimental results 

4.1 Mirror Reflectivity and Phase 

As described above, the ideal spectrographic calibration 
source spans more than an octave of spectral bandwidth; 
however, a cavity (which is necessary to filter the dense 
spectral modes) has limited spectral range over which its 
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modes overlap the evenly spaced comb modes. Let us first 
consider how much spectral bandwidth can be transmit- 
ted through a cavity. Figure [sj^a) depicts the experimental 
(green, solid) and theoretical (grey dashed) cavity-filtered 
spectrum normalized to the input spectrum while the cav- 
ity is locked at 980 nm. The experimental data are taken 
on a high-resolution optical spectrum analyzer with reso- 
lution bandwidth of 0.2 nm. The cavity effectively filters 
nearly 100 nm of the input spectrum before the cavity 
modes walk off the filtered comb modes. The fuzzy ap- 
pearance of the transmitted spectrum results from the op- 
tical spectrum analyzer sampling the discrete transmitted 
comb frequencies (and the absence of neighboring modes). 
The theoretical curve of lout/hn of each m}^ comb mode 
was generated by use of experimental parameters, mirror 
coatings of Fig.[2l and Eq. (T) Although the mirror reflec- 
tivity bandwidth appears to cover much of the spectrum 
from 800-1000 nm, the usable bandwidth from the cavity 
is limited by the variation of the phase shift of the mirrors, 
as discussed above. 

The overlap of the cavity and comb modes at the lock- 
ing wavelength of the cavity greatly affects the relevant 
bandwidth of the cavity filter, as shown in Figure [sjb). 
Here the power/mode of the comb spectra is shown, and 
the input curve is scaled by m = 20 to account for the fil- 
tered modes. The cavity parameters are initially optimized 
by use of the entire spectrum, and then different optical 
bandpass filters are placed in the path of the servo detec- 
tor. Optical spectra at a number of locking wavelengths 
are depicted. Note that the spectral width of the trans- 
mitted modes decreases as the cavity lock point is shifted 
away from 960 nm. This decrease results from the cavity- 
comb misalignment, or a shift in Av^c far away from the 
well aligned center modes. 

4.2 Cavity-Comb Offset 

If the the offset frequency fo of the laser is adjusted (by 
changing the dispersion in the laser cavity) such that the 
Fabry-Perot cavity length is an exact integral number of 
comb modes, but cavity and comb modes do not overlap 
at cJmid (ie., Avcc = /rep/2), it would seem that no light 
could be transmitted through the cavity. Without any cav- 
ity length adjustment, the comb modes are not in line with 
the cavity modes, and they are largely filtered. In prac- 
tice, however, the length of the cavity may be adjusted, 
very narrowly changing the cavity length away from an 
integral number of comb modes; in which case, the cavity 
length shift aligns the center frequency modes, but at the 
same time forces the cavity length to differ from an ex- 
act integral of comb modes. Figure |4] depicts theoretical 
curves adjusting Aucc near the center of reflectivity of the 
mirror coatings and plots the effect on the length adjust- 
ment of the cavity, total transmission through the center 
of the cavity, and the spectral bandwidth, at m = 20. To 
separate differing effects, the profiles assume a cavity in 
vacuum. 

As Aucc is adjusted, the cavity length has a corre- 
sponding linear adjustment as seen in Fig.ga). This ad- 
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Fig. 4. Calculation of the effect of Aucc of the comb relative to 
the cavity for the mirror design of Fig.|2] (a) As the offset shifts 
over the range of /rep(l GHz), the cavity length for optimum 
transmission and bandwidth vary linearly; at Aucc = /rep, the 
optimum alignment point shifts to the next peak in the length 
scan as represented by the discontinuities, (b) Less than 15 % 
variation in cavity transmission in a 50 nm bandwidth is seen 
as Ai^cc is varied. If mirrors of higher reflectivity are used, 
this variation increases, (c) Spectral width slightly increases as 
Ai^cc increases; however, note that the increase in bandwidth 
is at the expense of transmission. 

justment from the optimal cavity length causes a cor- 
responding decrease in the peak amplitude transmitted 
through the cavity. Figure is plotted using the ex- 
perimental mirror parameters. Higher mirror reflectivity 
would cause further reduction of the power through the 
cavity as the offset Aucc increases. For maximum trans- 
mission, the offset of the laser must be tailored to the 
cavity. While taking experimental data, Aucc is not opti- 
mized, but the cavity length is adjusted to compensate. 

It is perhaps counterintuitive that the spectral band- 
width (as defined by /out/^in > 0.5) passing through the 
cavity slightly increases as Aucc is increased [Fig.Qc)]. 
The broadening in the spectral width comes at the ex- 
pense of peak amplitude through the cavity in the case of 
vacuum. It should be noted that the spectral bandwidth 
may be increased further when accounting for the index of 
air in the cavity. The nearly-linear increase in index some- 
what balances the effect of the cavity mode's slip from the 
comb mode. 



4.3 Side-Mode Suppression 

Until now, the discussion of functional cavity bandwidth 
was centered on cavity transmission where the m^^ comb 
modes have >50 % transmission. For applications such as 
high-precision spectrographic calibration, the definition of 
usable bandwidth depends more critically on the suppres- 
sion of outlying frequency modes. An ideal calibration 
source resolves each m^^ comb mode. If other comb modes 
are not adequately suppressed, they can shift the spectral 
calibration causing systematic errors. 
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Fig. 5. (a) Heterodyne beat note of laser against filtered comb 
normalized to the m*^ mode. RF spectra shows both the cumu- 
lative RF beat (*) as well as the heterodyne (□) and conjugate 
heterodyne (o) beat notes, (b) Cycloid-like function of cumu- 
lative mode suppression as a function of RF frequency taken 
with the comb locked at different wavelengths. Red triangles 
[920 nm], black circles [960 nm] and blue boxes [1000 nm] repre- 
sent data taken in differing spectral regions. The red-dashed, 
black-dotted, and blue-dashed lines are theoretical simulations 
using experimental parameters. 



The standard method of spectral filtering has been 
demonstrated by measuring the RF spectrum of comb af- 
ter the filter cavity [26 . A fast photodiode detects the 
cumulative RF signal of all comb modes beating against 
each other. For a /rep = IGHz laser through a m = 10 filter 
cavity, the data and theory are shown in Fig.jsja) by the 
points marked by *. The first peak at 1 GHz results from 
the vector sum of each comb mode beating with its nearest 
neighbors; the second peak at 2 GHz from second nearest 
neighbors. The RF signal at 10 GHz derives from the sum 
of every 10*^ cavity mode. Let us examine the 1 GHz beat 
note. Taking into account the phase of the cavity as shown 
in Fig.[l] each individual 1 GHz beat note of a single comb 
mode against its neighbor will vary in phase from one to 
the next. The cumulative RF spectrum, each beat note 
a vector sum, may appear to have increased side-mode 
suppression. In the case where the comb and the cavity 
are not well aligned, the overestimate of suppression can 
be quite large. Figure Isfb) shows the cumulative suppres- 
sion measured with the cavity locked at 920 nm, where 
the comb and cavity are well matched, and at 960 nm and 
1000 nm where the cavity-comb alignment degrades. The 
theoretical curves (solid lines) take into account the vector 
sum of all beat notes. 

A more accurate measure of nearest neighbor suppres- 
sion may be achieved by heterodyning the output of the 
cavity against another laser. Rather than consider the sup- 
pression of all m — 1 filtered comb modes, here we are 
concerned only with the largest side modes. Due to the 



Airy function amplitude of the filter cavity, the modes 
which are least suppressed are those directly neighboring 
the m^^ comb modes. Figure [I] defines the left- and right- 
nearest- neighbor frequency modes (nni and nur) at the 
center frequency of the spectrum. Because nni and nur 
modes have the least suppression of the filtered modes 
(unless the cavity and comb are completely misaligned), a 
more meaningful definition of the usable bandwidth may 
be defined using the suppression of these modes. Using the 
amplitude from Eq. [l] for the filtered modes with respect 
to the m^^ modes, the maximum suppression S in dB of 
the k^^ sideband at reflectivity R is: 



S = 10 log 



2 cos 



\ m/rep J 



R^l 



R^ -2 cos 



27r 



\ mf rep m J 



R^l 



(2) 



where = ±1 refers to the nearest neighbor modes. 

The □ beat notes in Fig.jsja) denote the RF frequencies 
resulting from a ~ 4 mW cw diode laser at 960 nm beating 
against the m^^ comb tooth. By measuring the suppres- 
sion of the right and left □ beat notes against the peak, a 
more physical measurement of the suppression of nearest 
neighbors (nni and nur) may be obtained. The conjugate 
beat note of the laser against the m + 1 mode is marked 
with o. To use the filtered comb as calibration source, it is 
necessary to know the suppression of intermediate modes 
as a function of wavelength. Because the heterodyne laser 
is not tunable across the entire bandwidth, the comb is 
locked by use of different frequency bands of the output 
cavity spectrum. Figure [6] shows the theoretical and ex- 
perimental curves at m = 10 (solid and filled circles) and 
m = 20 (solid and filled diamonds) for the suppression of 
nni and nn^. The simulated curves take into account the 
full theoretical model where the length of the cavity is ad- 
justed as a free parameter. Both m = 10 and m = 20 show 
similar dependence upon locking wavelength; however, the 
maximum suppression becomes greater as the number of 
intermediate modes to be filtered is decreased. In a high- 
resolution spectrograph where side-mode suppression can 
cause an apparent shift in calibration, the limiting fac- 
tor for the usable spectral bandwidth of the cavity results 
from both the suppression of the side modes as well as 
from the asymmetry of the suppressed modes. 



4.4 Filtering an Octave of Comb Bandwidth 

Use of a frequency comb for astronomical spectrograph 
calibration requires (1) a large spectral bandwidth, (2) 
frequency modes spaced by 10-30 GHz, and (3) high sup- 
pression (>27dB) of neighboring modes for velocity cal- 
ibrations accurate to cm/s [1^. These combined require- 
ments are unphysical for a single filter cavity with our 
current coatings; however, multiple cavities in parallel al- 
low coverage over the entire spectrum [2]. As depicted in 
Fig.[7| the comb spectrum may be split into functional 
spectral bandwidths by use of dichroic mirrors. Each cav- 
ity (cl-c4) filters the comb to m/rep, and the filtered combs 
are recombined. The final mirror, Mfn, may be designed 
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Fig. 6. Data are plotted by measuring the suppression of the 
nni and nur at each of the locking bandwidths. Open and filled 
diamonds are the suppression of nni and nur at m = 20, and 
open and filled boxes are for m = 10. The solid lines represent 
the theoretical curves for right (solid line) and left (dashed 
line) nearest neighbors. 



to spectrally flatten the comb output. In place of Mfii, a 
grating/spatial light modulator system could be used. Al- 
though the filtered frequencies of each cavity are spaced 
by m/rep, the cavities may have a frequency offset. For 
example, cavity cl may filter the m modes over the blue 
wavelengths while cavity c2 could pass the m + 5 modes 
{i.e. both cavities have a filter number m, but may pick 
out different comb teeth). This offset, in which the m^^ 
teeth are filtered, can then be used to reference one cavity 
to the next. In the above example, the beat note between 
the overlapping spectral region of cl and c2 is 5/rep. 




Fig. 7. Because a single cavity with mirror coatings of Fig.|2| 
cannot filter a wide spectrum, the comb is spectrally divided 
among a number of filter cavities. Cavities (cl-c4) process the 
spectrum in parallel, each cavity filtering a portion of the spec- 
trum to the necessary frequency-mode spacing. The length of 
individual cavities is stabilized by a PZT, and the offset be- 
tween the cavities (an integral of /rep) is monitored with de- 
tectors (dl2, d23, d34). 



The offset of each cavity with respect to the next may 
be designed directly into the mirror coatings, or more 
practically, a slight change in cavity length may be used at 
the expense of cavity transmission bandwidth. For more 
control, a gas of known index may be used in combination 
with a small change in the length of the cavity. 



5 Spectrograph Calibration 

Until now, we have considered only the frequency filtering 
of comb modes to spectrographically resolvable frequency 
spacing while maintaining suppression of the neighboring 
frequency modes. For this comb to be a useful calibration 
source, we must also address precision. Although the con- 
stant spacing of the comb lines has been demonstrated, 
many factors can create an apparent shift of the center of 
gravity (COG) of the comb lines. 

Figure [S] shows an image of 35 comb lines that have 
been filtered by a 20 GHz cavity, dispersed by a high- 
resolution grating, and imaged onto a charged-coupled de- 
vice (CCD) camera. The graph plots a subset of CCD pix- 
els as a function of intensity and shows nine of these comb 
lines. For a standard spectrographic calibration, the data 
would be fit with nine individual Gaussian functions, each 
with amplitude, width, position, and offset. The strength 
of using a comb to calibrate the spectrograph lies in the 
ability to leverage the equally spaced frequency compo- 
nents. Small errors that would typically cause an apparent 
shift in the center of gravity of one calibration point are 
now averaged out. In place of N Gaussians, each with four 
free parameters, the entire spectrum may be fit with a sum 
over N Gaussians, equally spaced along x. Here the sum 
aExp[— c(x — bN — d)^] + e uses only five free param- 
eters for the entire comb: peak a, separation 6, width c, 
uniform shift (i, and noise floor e. Because we assume con- 
stant spacing, the fit is less sensitive to parameters such 
as asymmetric suppression of nearest-neighbor modes. 
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Fig. 8. (a) CCD image of a 1 GHz Mode-locked laser filtered 
to 20 GHz by a cavity and then spectrally dispersed with a 
grating, (b) Intensity plot of subsection. 

Although calibration in this manner does reduce some 
systematic errors, the properties of the comb itself must be 
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considered. In the above simulations, the comb hnes were 
assumed to be infinitely narrow. Under this premise, the 
filter cavity acts on the amplitude of a delta function, and 
no COG shift is possible. In a physical comb system, the 
frequency reference, which is used to stabilize the comb, 
contributes to a linewidth of the comb teeth. When the 
comb is locked, the frequency modes of the laser may be 
modeled as a Gaussian with a finite full width at half max, 
z^L- If the cavity and comb modes are perfectly overlapped, 
the comb spectral distribution may slightly change the 
shape of the comb mode, but does not shift the COG. 
As the cavity and comb slip off one another, the cavity 
asymmetrically reshapes the comb teeth, thereby giving 
the appearance of a frequency shift. The larger the ratio of 
comb tooth width to the cavity width, the greater the shift 
in COG. In the limit where the width of the comb tooth 
exceeds that of the cavity, the cavity function completely 
reshapes the comb mode, and there is a 1:1 ratio between 
the COG shift of the comb line and Aucc- This systematic 
error would be carried through the spectrograph, creating 
an apparent shift in calibration. Figure [9] plots the COG 
shift as a function of the offset between cavity and comb 
for m = 20 (a) and m = 10 (b). These curves are plotted 
for a mirror refiectivity of R = 0.99. On the upper axis, 
Eq.[2j is used to denote suppression of the nearest neighbor 
with respect to the m^^ comb mode. The combination 
of necessary suppression and maximum tolerable shift in 
COG determine the allowed linewidth of the comb modes. 
In a cavity of R = 0.99 and m = 10, which results in a 
cavity full- width at half-max of 32 MHz, a comb linewidth 
of Aui^ = 1 MHz has less than 20 kHz shift of the COG 
while suppression is better than 30 dB. With this level of 
stabilization, spectroscopic calibrations on the level of a 
few centimeters per second are possible. For cavity mirrors 
with refiectivity greater than 99%, the shift in COG is 
greater than that plotted in Fig.|9] 

Depending on the level of precision necessary for spec- 
trograph calibration, various frequency reference may be 
used to stabilize the comb. A practical and relatively in- 
expensive option is a low-noise quartz oscillator (typically 
5 or 10 MHz) that has its frequency steered on interme- 
diate time scales to a compact microwave atomic stan- 
dard (such as Rb). On longer time scales the accuracy of 
the frequency reference can be guided by signals from the 
constellation of satellites that form the global position- 
ing system (GPS). Such frequency references are commer- 
cially available and often referred to as GPS-Disciplined 
oscillators (GPSDO) [27'28'. The fractional frequency un- 
certainty available from a GPSDO can be near or below 
the level of 1 x 10~^^ for averaging times greater that 1 
second. The uncertainty decreases approximately propor- 
tional to the averaging time on scales greater than ~ 1 
day, making feasible high-accuracy measurements even if 
they are separated by days or years. 

It is important not to confuse the frequency uncer- 
tainty of the GPSDO with the linewidth it would pro- 
vide for the frequency comb elements, as we might er- 
roneously assume that an uncertainty of 1 x 10~^^ pro- 
vides a linewidth of a few kilohertz for the optical comb 



elements. In fact, the determination of the exact comb 
linewidth can involve many factors including, (1) the fre- 
quency noise on the free-running femtosecond laser fre- 
quency comb, (2) the short-term (< Is) phase-noise of 
the quartz oscillator in the GPSDO, (3) the residual phase 
noise on a synthesizer used in converting from the quartz 
frequency to the repetition rate of the femtosecond laser, 
and (4) the electronic servo system used in controlling 
the femtosecond laser comb relative to the GPSDO. The 
various contributions of these elements to the observed 
linewidth of the comb elements can be calculated, keep- 
ing in mind that the noise on the quartz crystal and mi- 
crowave electronics is multiplied by the ratio of the optical 
and radio/microwave frequencies {e.g., something between 
5 X 10^ and 5 x lO*"). This large multiplication factor can 
transform the apparently low noise of a 10 MHz oscillator 
(an the intermediate microwave synthesizer) into an opti- 
cal linewidth that is much greater than we would naively 
anticipate. Nonetheless, with reasonable care and read- 
ily available components, a linewidth of the comb modes 
that is less than 1 MHz in the visible portion of the spec- 
trum is readily attained. Improved radio and microwave 
frequency electronics could further reduce this linewidth. 
As just noted, the linewidth scales with frequency, so a dif- 
ference in linewidth by a factor of ~ 2 could be expected 
in the two extremes of an octave-spanning spectrum. If 
a still narrower linewidth for the optical comb modes is 
required, at the expense of simplicity we could employ a 



(a) 




29 




NN suppression [dB] 
28 27 26 


25 24 


[kHz] 








' Avl=10MHz 


— 1 -1- 
















G shift 








Avl=1 MHz 


















O 














^ 0.01 








, Av =0.1 MHz 
















10 


20 30 
Av^^ [MHz] 


40 


(b) 

1000 
X 100 
^ 10 

^ 1 

o 

O 1 




35 




NN suppression [dB] 
34 33 32 


31 








— Avl=10MHz 


1 






















■ Avl=1 MHz 










u 


^ Av, =0. 1 MHz 

^ 1 1 ^ . 1 










5 


10 15 
Av,, [MHz] 


20 



Fig. 9. Apparent COG shift of a filtered comb line as a func- 
tion of the offset between the cavity and a comb mode. Sup- 
pression of the nearest neighbor mode is plotted on the upper 
axis. Reflectivity of = 0.99 is assumed, and the shift is plot- 
ted for Avi^ of 0.1 MHz, 1 MHz, and 10 MHz. Filter numbers 
of m = 20 and m = 10 are shown in (a) and (b), respectively. 
For a /rep = 1 GHz comb and R — 0.99, this corresponds to a 
cavity linewidth of 64 MHz at m = 20 and 32 MHz at m = 10. 
At ^800 nm, a shift of 12.5kHz is equivalent to a radial velocity 
of 1 cm/s. 
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stabilized laser as a short time (< 1 s) reference in place 
of the quartz oscillator. In such a case, optical linewidths 
from 100 kHz down to ^1 Hz can be achieved [29 . 



6 Conclusion 

Frequency combs are a promising avenue for the calibra- 
tion of astronomical spectrographs. Current comb technol- 
ogy requires filtering of comb modes to produce adequate 
mode separation that can be resolved by standard spec- 
trographs. The parameters relevant to cavity filtering have 
been detailed with a demonstrated agreement between ex- 
periment and theory, and factors that impact COG line 
shifts are discussed. Beyond the compelling benefits of fre- 
quency combs as a calibration source for observational as- 
tronomy and cosmology, a laser frequency comb calibra- 
tion source could positively impact several additional ar- 
eas that require robust frequency sources with 10-20 GHz 
mode spacing. Examples include direct spectroscopy with 
frequency combs, remote sensing, optical and microwave 
waveform synthesis, high-speed coherent communications, 
and optical clock development. 
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design. Finally, we thank Nathan Newbury and Elizabeth 
Donley for their thoughtful comments on this manuscript. 
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